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摘要 : 优化 台 址 风 环境 对 提高 大 口径 、 高 精度 射电 望远镜 在 中 高 频段 的 有 效 观 测 时 长 具 
有 重要 意义 。 通 过 风 障 调控 风 场 可 以 有 效 减 小 风 和 荷载 对 望远镜 的 影响 。 风 障 的 布局 设计 除了 
与 风 障 高 度 、 孔 阶 率 等 参量 有 关外 ， 还 要 综合 考虑 台 址 地 形 对 挡 风 效率 的 影响 。 以 奇 台 射电 
望远镜 (QiTai radio Telescope, QTT) 台 址 为 研究 对 象 ， 开 展 不 同 地 形 对 风 场 调控 影响 的 数值 
模拟 研究 。 仿 真 结果 表明 ， 无 风 障 时 望远镜 区 的 风速 大 小 主要 与 上 下 游 边界 地 形 的 最 高 海拔 
有 关 ; 布设 风 障 时 ， 若 望远镜 区 域 为 水 平地 形 ， 风 障 遮 项 区 降低 风速 的 大 小 与 上 游 入 流风 攻 
角 以 及 下 游 地 形 的 海拔 高 度 有 关 ; 台 址 实际 斜坡 地 形 增加 了 流 场 的 复杂 程度 ， 在 进行 风 障 风 
场 仿真 时 ， 上 游 边 界 要 尽量 延伸 到 相对 较 高 山体 的 外 围 。 不 同 边界 地 形 的 仿真 结果 表明 ， 望 
远 镜 区 风速 折 减 效率 差 值 最 大 达到 6%。 该 研究 可 以 为 望远镜 台 址 风 场 调控 仿真 中 的 地 形 建 


T 模 提 供 可 靠 的 参考 。 
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BR. WLU ABS 


Bog, XPS TD RE REOR ies), RAT PEP PSP LU, dE 
向 精度 不 可 避免 地 受到 台 址 风 扰 影响 “” 。 对 于 大 口径 、 高 指向 精度 射电 望远镜 而 言 ， 风 扰 影响 天 线 
© Ji rap E BE RS [s AN AL), RRE, HE tl eb ee e CT Be E JEA Ui AY BA A e R BE 
= 抗 风 需 求 '” , ER SET EE, AUP oI 7752 EEE, AAD A EIT 


£z 防风 抑 侍 '”; AUR ELS AR SAS, EB SETTE e E DARLES E P DLS ^7 156 ole D EE 00 £5 e cU; 
Q URK AAA BEER RUZI 。 风 障 是 一 种 具有 一 定 开 孔 率 或 者 完全 不 透风 的 障碍 物 ， 对 
风 可 以 产生 一 定 的 阻碍 作用 ， 从 而 改变 风向 并 消耗 风能 ， 使 得 障 后 的 风速 衰减 。 利 用 风 障 调控 风 场 一 

定 程度 上 可 以 改善 台 址 的 风 环 境 。 

目前 ， 针 对 风 障 调控 风 场 的 研究 主要 集中 在 风 障 孔隙 率 、 高 度 、 面 形 等 参量 对 挡 风 效率 的 影响 方 

面 ': 一 ”。 由 于 布设 风 障 领域 的 地 形 通常 相对 比较 平坦 ， 研 究 中 对 地 形 的 关注 度 比较 小 。 而 大 口径 、 

高 精度 射电 望远镜 台 址 大 多 位 于 山谷 或 倪 地 内 中 ， 地 形 相对 复杂 ， 风 的 流 场 结构 也 与 其 他 领域 的 情况 

有 所 不 同 。 文 [18] 研 究 了 不 同 山 丘 地 形 的 风 场 绕 流 分 布 ， 结 果 表 明 高 宽 比 较 大 的 山 丘 背风 面 绕 流 稳定 性 


比较 差 ， 易 发 生 流动 分 离 ; 山 丘 越 高 ， 顶 部 加 速 比 越 大 ， 下 游 的 尾 流 恢复 距离 


工 长 。 文 [19] 利 用 风 洞 试 


验 和 数值 模拟 方法 研究 了 某 山地 地 形 的 风 场 特性 ， 指 出 山地 风速 、 风 向 变化 剧烈 ， 山 峭 位 置 的 风速 比 在 
1.6-2.8 之 间 ， 山 谷 位 置 的 风速 比 在 0. 12~1.6 之 间 。 文 [20] 基 于 缩 比 模型 开展 了 山地 环境 的 风 洞 实 
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验 ， 研 究 山 地 背风 面 的 尾 流 分 布 特征 ， 人 研究 结果 表明 ， 山 体 下 游 的 流 场 结构 与 背风 坡 面 角度 有 关 。 

为 研究 山地 环境 中 地 形 对 大 口径 射电 望远镜 台 址 风 场 调控 的 影响 ， 本 文 以 奇 台 射电 望远镜 台 址 地 
形 为 研究 对 象 ， 利 用 数值 模拟 方法 开展 不 同 地 形 对 风 场 调控 影响 的 研究 。 以 天 线 位 置 为 中 心 划分 核心 
区 域 ， 核 心 区 域 的 两 端 分 别 设置 风 障 ， 在 核心 区 域 的 南北 方向 分 别 设置 3 组 不 同 地 形 的 边界 。 为 了 比 
较 台 址 斜坡 地 形 对 风 场 调控 的 影响 ， 构 建 台 址 水 平地 形 对 照 组 ， 通 过 开展 不 同 地 形 工 况 的 风 场 仿真 研 
究 掌握 地 形 对 风 场 调控 的 影响 关系 ， 为 风 障 优化 布局 提供 参考 。 


1 数值 模拟 数学 模型 


1.1 庙 流 控制 方程 

本 文采 用 雷诺 平均 (Reynolds Average Navier-Stokes, RANS) 数值 模拟 方法 进行 仿真 。 对 于 野外 小 
区 域 风 场 仿真 ， 可 以 将 空气 假设 为 低速 、 不 可 不 的 粘性 牛顿 流体 。 风 的 流动 遵循 质量 守恒 和 动量 守恒 
定律 ,控制 方程 为 


op ð 
十 .)=0 1 
ar tag, (Om)=0, (1) 


9(puj) 9 dp ð du, Ou, 2 Ou; ð 一 一 
2 quisi soe ROPA M T ES ES 
d jg PR ae ae ax, x, 3 ax] dx, eS (2) 


其 中 , p 为 空气 密度 ; 1 为 时 间 ; x 为 沿 i 方 向 的 坐标 ; w 为 i 方向 的 速度 ; p 为 压力 ; u 为 流体 粘度 ; 
© 9,79 Kronecker PAR; — puu! 为 雷诺 应 力 ，5; 为 i 项 动量 源 。 

i HET EASE ALA, PES Aa LET RF, MARASKA TA 
求解 精度 更 高 的 可 实现 k-e 模型 ， 其 滑动 能 上 和 耗 散 率 e 的 输 运 方程 为 


ðk ð M, \ ðk 
E +G. +G, -pe - Y, , 3 
Po eus) ud k b TPE M (3) 
aa a A — +0. 5 c.6 (4) 
- — 
p ðt Ox; H T, Ox; p T p ^q i: TE le k 3e 7b 


JUP, o, 为 滑动 能 上 MUTE DIRE, o, 为 耗 散 率 e MINGIS C, 为 由 平均 速度 梯度 引起 的 
滑动 能 产生 项 ;G 为 由 浮力 引起 的 满 动能 产生 项 ，Y, 为 可 压缩 淇 流 脉动 对 总 耗 散 率 的 影响 ; 满 流 炸 
O HERM = pC, C, =max[0.43, T| ,= Sh/e ;模型 中 的 一 些 党 数量 由 经 验 值得 到 ,C=1. 4 
C,=1.92, og,=1.0, o,=1.2, 
1.2 多 了 筷 跳 路 数学 模型 

本 文 的 风 障 具有 一 定 的 孔隙 ， 可 以 看 作 多 孔 介 质 ， 仿 真 建 模 可 以 采用 多 孔 跳跃 模型 。 模 型 构建 原 
理 是 在 动量 方程 中 增加 动量 源 项 。 该 源 项 由 黏 性 损失 项 和 惯性 损失 项 两 部 分 组 成 ， 方 程 为 


1 
s, =- (Eu «e T olus). (5) 
a 2 


其 中 ,a 为 渗透 率 ; C. 为 惯性 阻力 系数 ;u 为 流体 速度 。 仿 真 中 改变 渗透 率 a 对 仿真 结果 的 影响 极 
小 ， 该 参量 可 以 设置 为 常数 。 惯 性 阻力 系数 C, 可 以 采用 经 验 公式 
_K _1-B 
C = 人 ， (6) der (7) 


HL ， 其 中 ， 开 为 压力 损失 系数 ; B 为 风 障 的 孔隙 率 ; An 为 风 障 的 厚度 。 本 文中 风 障 孔隙 率 设 
置 为 33% ， 风 障 厚度 设置 为 0.001 6 m, 
1.3 大气 边界 层 边 界 条 件 

大 气 边界 层 内 ， 风 速 随 高 度 变化 的 规律 可 以 采用 对 数 律 或 指数 律 数学 公式 描述 。 计 算 域 人 口 风 速 
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采用 更 适用 于 工程 计算 的 指数 律 公 式 表 达 : 
zo 
V(z) = A : (8) 
Zo 


台 址 地 貌 类 型 属于 B 类 粗糙 度 ，a 指数 取 0. 15; z 为 海拔 高 度 ; z CIO m; VW 取 15 m/s, 
计算 域 入 口 满 流 强度 采用 日 本 建筑 规范 中 卫 类 地 貌 的 满 流 强度 经 验 公式 


0.23, z € 5 
S -0.2 9 
Ae oaf) TI |: E 
z 
g 
z, 取 350m, IIRA Dab aE k Rin D AE BCR e 按照 
3 ; B 
k => Volo, (10) e = 0.09 P (11) 


设置 ， 其 中 ，Z 为 清流 积分 尺度 ， 取 0. 07 倍 的 计算 域 项 高 。 
2. WU ESAME 


2.1 奇 台 射电 望远镜 台 址 计算 域 模型 构建 

奇 台 射电 望远镜 台 址 位 于 一 处 四 面 环 山 的 矩形 
盆地 内 ， 地 理 坐 标 为 43?36'4. 03"N, 89^40'56. 99"E , 
海拔 约 为 1760m， 台 址 地 形 如 图 1。 台 址 园区 内 
的 来 风 主要 在 南北 方向 ， 东 西方 向 来 风 较 少 [31。 
因此 ， 本 文 主要 围绕 南北 方向 的 地 形 对 台 址 风声 
调控 产生 的 影响 开展 仿真 研究 。 

风 场 调控 数值 模拟 可 以 采用 二 维 或 三 维 仿 
真 ， 针 对 单一 变量 的 研究 采用 二 维 仿 真 优势 更 为 
明显 22。 本 文采 用 二 维 仿真 ， 沿 天 线 所 
在 点 的 南北 方向 提取 台 址 地 形 高 程 数 据 ， 如 图 
2。 以 天 线 位 置 为 中 心 ， 将 水 平方 向 -600~ 600 m 
的 区 域 划分 为 核心 区 域 ， 以 核心 区 域 为 中 心 分 别 
向 北向 南 延 伸 。 向 北 延伸 的 端点 编号 为 1, 2, 3; 
向 南 延 伸 的 端点 编号 为 a, b, c。 将 2 组 边界 端 
点 分 别 进行 组 合 ， 共 产生 9 组 地 形 模型 (1a, 1b, 


ak 
f= 
PE 
aa 
LJ 


le, =, 3a, 3b, 3c) 。 以 地 形 模型 为 基础 构建 台 图 1 Foe Bt SHH 
址 风 场 仿真 计算 域 模型 ， 计算 域 模型 如 图 3。 在 Fig. 1 The terrain of the QTT site 


有 风 障 的 工 况 中 ,分别 在 水 平方 向 -600 m 和 600 m 
的 位 置 设置 风 障 ， 风 障 高 度 设置 为 120 m。 计 算 域 模型 顶 高 设置 为 3000m。 为 了 保障 所 有 仿真 计算 域 
入 口 风 速 相同 ， 台 址 实际 地 形 以 外 做 了 水 平缓 冲 和 地 形 缓冲 。 水 平缓 冲 相 对 天 线 位 置 的 高 度 和 长 度 在 
所 有 仿真 计算 域 模型 中 相同 。 来 风 条 件 设置 为 南方 向 来 风 和 北方 向 来 风 。 
2.2 ” 台 址 实际 斜坡 地 形 仿真 结 

各 工 况 仿真 结果 提取 相同 位 置 的 风速 数据 。 竖 直方 向 上 ， 提 取 相 对 地 面 y,= 60 m(0.5 倍 的 风 障 
高 度 ) 96 m(0.8 倍 的 风 障 高 度 ) 高 度 的 数据 ; 水 平方 向 上 ， 从 -600~600 m, fifi 50 m 提取 一 个 数 
据点 。 风 障 的 风速 折 减 系数 采用 


| x 100% (12) 


HSE, 其 中 ,A 为 风速 折 减 系数 ; V. 为 无 风 障 工 况 风速 ;WV, 为 有 风 障 工 况 风速 。 
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图 2 台 址 地 形 高 程 数据 


Fig.2 Terrain elevation data of the site 


台 址 实际 斜坡 地 形 顺 风向 的 风速 折 减 系数 分 7 
布 如 图 4。 由 图 4 可 以 看 出 ， 北 方向 来 风 ， 顺 风 Outlet 
向 风速 折 减 系数 从 -600 m 位 置 的 风 障 后 先 增 大 ， ”|/ N 
然后 减 小 在 600 m 位 置 的 风 障 前 再 次 增 大 ; 南 Buffer ofterrain 
方向 来 风 ， 顺 风向 风速 折 减 系数 从 600 m 位 置 的 E J NE N 
风 障 后 近似 持续 增加 。 上 游 边界 同 起 点 的 3 组 地 x 
形 工 况 ， 风 速 折 减 系 数 曲 线 变 化 更 为 接近 ， 如 图 
4(a) 中 上 游 起 点 为 “1” 的 3 组 数据 。 由 于 下 游 
边界 不 同 ， 曲 线 之 间 略 有 差 值 。 由 此 说 明 上 游 边 
界 地 形 对 风速 折 减 系数 的 影响 大 于 下 游 边 界 地 形 的 影响 。 天 线 位 置 (水 平 0 点 ) 的 气流 相 比 于 风 障 位 
置 更 加 稳定 ， 提 取 天 线 位 置 相对 地 面 y,=96 m 高 度 处 在 有 风 障 工 况 和 无 风 障 工 况 的 风速 以 及 各 工 况 的 
风速 折 减 系数 如 表 1。 为 了 精确 比较 不 同上 下 游 边 界 地 形 对 风 场 调控 的 影响 ， 制 作 地 形 工 况 1c, 2c, 
3c, 3b, 3a 在 有 风 障 工 况 和 无 风 障 工 况 的 顺风 向 风 场 流 线 图 ， 如 图 5 和 图 6。 

结合 图 4、 图 5 以 及 表 1 分 析 北 方向 来 风 边 界 地 形 对 风 场 调控 的 影响 。 上 游 边界 不 同 、 下 游 边 界 
相同 地 形 的 比较 (地 形 1c, 2c 和 3c) ， 在 无 风 障 工 况 中 上 游 边界 内 最 高 点 的 海拔 高 度 越 高 ， 到 达 天 线 
位 置 的 风速 越 小 。le 地 形 上 游 最 高 海拔 212 m， 天 线 位 置 的 风速 为 23. 63 m/s; 3e 地 形 上 游 最 高 海拔 
28 m， 天 线 位 置 的 风速 为 24. 02 m/s。 在 有 风 障 工 况 中 上 游 边界 内 最 高 点 的 海拔 高 度 越 高 ， 经 风 障 调 
控 到 达 天 线 位 置 的 风速 反而 越 大 ，lec 地 形 天 线 位 置 的 风速 为 13. 62 m/s, 3c 地 形 天 线 位 置 的 风速 为 
11.75 m/s。 结 合 风 场 流 线 图 分 析 认 为 ， 无 风 障 工 况 中 ， 上 游 山 体 越 高 ， 对 风 的 遮挡 能 力 越 强 ; 而 有 
风 障 工 况 中 ， 风 障 在 天 线 位 置 削弱 风速 的 大 小 与 进入 核心 区 域 的 风 攻 角 呈 正 相关 ( 无 风 障 时 -600 m 位 
置 的 风 攻 角 ) ，le, 2e 和 3c 进入 核心 区 域 的 风 攻 角 分 别 为 2.28”，2.73 和 3. 09 。 下 游 边 界 不 同 、 上 
游 边界 相同 地 形 的 比较 (地 形 3c, 3b 和 3a) ， 无 风 障 工 况 中 ， 天 线 位 置 的 风速 同样 与 下 游 边界 内 最 高 
点 的 海拔 高 度 呈 负 相 关 。 有 风 障 工 况 中 ，3e 地 形 下 游 边界 的 地 形 遮 挡 最 小 ， 但 由 于 风 出 口 处 发 生 了 
流动 分 离 ， 产 生 回 流 ， 对 风 的 流动 形成 阻碍 (如 图 5“3c-N-F”) ， 使 得 天 线 位 置 的 风速 最 小 。3b 地 形 
在 风 出 口 处 未 发 生 流动 分 离 ， 且 地 形 阻碍 比 3a 小 ， 在 天 线 位 置 的 风速 最 大 。 

结合 图 4、 图 6 以 及 表 1 分 析 南 方向 来 风 边 界 地 形 对 风 场 调控 的 影响 。 无 风 障 工 况 中 ， 天 线 位 置 
的 风速 均 与 上 下 游 边界 内 最 高 点 的 海拔 高 度 呈 负 相 关 。 有 风 障 工 况 上 游 边界 不 同 、 下 游 边界 相同 地 形 
的 比较 (地 形 3c, 3b 和 3a) ，3 组 地 形 在 600 m 位 置 的 和 流风 攻 角 相近 ,但 由 于 和 斜坡 地 形 和 风 障 的 共 
同 作用 ， 天 线 位 置 发 生 了 流动 分 离 ， 被 回流 旋涡 覆盖 。 从 风 场 流 线 图 中 明显 可 以 看 出 ， 上 游 风速 越 
大 ， 旋 涡 的 范围 越 大 ， 而 天 线 位 置 的 风速 与 入 流风 速成 正比 。 图 4(c) 中 上 游 起 点 为 “ec” 的 3 组 数据 
在 0 点 下 游 风 速 折 减 系数 开始 下 降 ， 而 图 4(d) 中 风速 折 减 系数 则 继续 增强 ， 从 风 场 流 线 图 中 可 以 看 
出 (图 6 中， 有 风 障 工 况 3c, 2e, lc 地形)， 旋 涡 中 心 更 接近 地 面 96 m 高 度 的 位 置 ， 风 速 相对 更 小 。 
而 地 面 60 m 高 度 距 离 旋 涡 中心 更 远 ， 风 速 相 对 较 大 ， 且 从 0 点 开始 ， 下 游 回 流速 度 明 显 增强 ， 因 此 
风速 折 减 系数 降低 。 有 风 障 工 况 中 ， 上 游 边界 相 同 、 下 游 边界 不 同 地 形 的 比较 (地 形 3c, 2c Mle), 3 
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图 3 仿真 计算 域 


Fig. 3 Computational domain of simulation 
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组 地 形 仿真 结果 天 线 位 置 的 风速 差 值 很 小 ,但 同样 受 回流 旋涡 位 置 的 影响 。 回 流 旋涡 形成 的 位 置 和 范 


围 受 风 障 自身 参量 、 斜 坡地 形 、 下 游 边界 地 形 契合 作用 影响 。 


一 一 3a 一 四 一 3b —xk— 3c 


s 
A 
-600  -400  -200 0 200 400 600 -600  -400  -200 0 200 400 600 
Distance from the antenna / m Distance from the antenna / m 
(a) (b) 
100 F 100 F 
90r 90 F 
80 F 80 F 
= 70 上 = 7+ 
< < 
60r 60 上 
Wind from south -— Wind from south -< 一 
50r y,760m 50r y,=96 m 
40 1 1 1 1 1 1 1 40LL_ LIL LUL | 
-600  -400  -200 0 200 400 600 -600  -400  -200 0 200 400 600 
Distance from the antenna / m Distance from the antenna / m 
(c) (d) 
图 4 ” 台 址 斜坡 地 形 顺 风向 风速 折 减 系数 分 布 
Fig.4 Distribution of reducing wind speed coefficients downwind on the slope terrain of the site 
R1 台 址 斜坡 地 形 天 线 位 置 风 速 折 减 系数 
Table 1 Ranking of reducing wind speed coefficients for antenna location on slope terrain of site 
Wind from north Wind form south 
Wind speed Wind speed Reducing Wind speed Wind speed Reducing 
Terrain without wind with wind wind speed Terrain without wind with wind wind speed 
number barrier barrier coefficient number barrier barrier coefficient 
/ (m/s) /(m/s) / (96) / (m/s) / (m/s) / (96) 
3c 24. 02 11. 75 51.08 le 23. 65 1.67 92. 94 
3b 23. 80 11. 89 50. 03 3c 23.98 1. 72 92. 82 
3a 23.61 11. 88 49. 68 2c 23.90 1.77 92. 61 
2c 23. 84 12.19 48. 85 1b 23.27 2.93 87. 42 
2b 23. 60 12. 32 47. 76 3b 23. 67 3. 09 86. 95 
2a 23. 40 12. 33 47.31 2b 23.54 3.11 86. 79 
le 23. 63 13. 62 42. 37 la 22. 24 3.51 84. 22 
1b 23.38 13. 75 41. 18 3a 22. 66 3.70 83. 65 
la 23.17 13. 75 40. 65 2a 22.59 3. 74 83. 46 
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图 5 北方 向 来 风 台 址 实际 斜坡 地 形 顺 风向 风 场 流 线 


Fig.5 Downwind wind flow streamline on the slope terrain of the site with northward wind 
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图 6 南方 向 来 风 台 址 实际 斜坡 地 形 顺 风向 风 场 流 线 


Fig.6 Downwind wind flow streamline on the slope terrain of the site with southward wind 
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2.3 台 址 水 平地 形 构建 及 仿真 结果 分 析 

核心 区 域 属 于 斜坡 地 形 ， 处 于 迎风 坡 与 背风 坡 时 的 风 场 结构 明显 不 同 ， 风 障 调控 风 场 的 效果 差异 
也 很 大 。 由 于 和 斜坡 地 形 的 影响 ， 图 4 中 ， 北 方向 来 风 ， 南 边 边界 属于 风 场 下 游 ，3 组 下 游 边界 地 形 最 
高 点 海拔 落差 相对 较 小 ,但 对 风速 折 减 系数 影响 的 差异 明显 。 而 南方 向 来 风 ， 北 边 边界 3 组 地 形 最 高 
点 的 海拔 落差 相对 较 大 ， 但 对 风速 折 减 系数 影响 的 差异 反而 相对 不 明显 。 为 了 比较 斜坡 地 形 对 风 场 流 


动 的 主导 影响 以 及 明确 不 同 边界 地 形 对 核心 区 域 风 场 调 控 的 影响 关系 ， 将 原 有 和 斜坡 地 形 构 建 为 相对 水 
平 的 地 形 进行 风 场 仿真 研究 ， 如 图 7。 仿真 计算 域 模型 构建 与 2. 1 节 相 同 。 
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图 7 台 址 水 平地 形 构建 


Fig.7 Constructing horizontal terrain of the site 


台 址 水 平地 形 顺 风向 的 风速 折 减 系数 分 布 如 图 8。 由 图 8 可 以 看 出 ， 台 址 斜坡 地 形 “ 变 ”水 平 后 ， 
北方 向 来 风 和 南方 向 来 风 的 风速 折 减 系数 变化 趋势 相同 ， 均 为 先 增 大 再 减 小 ， 然 后 增 大 。 提 取 天 线 位 
置 相 对 地 面 y,=96 m 高 度 处 在 有 风 障 工 况 和 无 风 障 工 况 的 风速 以 及 各 工 况 的 风速 折 减 系数 如 表 2。 制 
作 地 形 工 况 1e, 2c, 3c, 3b 和 3a 在 有 风 障 工 况 和 无 风 障 工 况 的 顺风 向 风 场 流 线 图 ， 如 图 9 和 图 10。 

根据 图 8 和 表 2， 无 风 障 工 况 中 ， 天 线 位 置 的 风速 均 与 上 下 游 边界 地 形 最 高 点 的 海拔 有 关 ， 上 下 
游 边界 地 形 最 高 点 海拔 越 高 ， 天 线 位 置 的 风速 越 小 。 结 合 图 8、 图 9 以 及 表 2 分 析 北 方向 来 风 边界 地 
形 对 风 场 调控 的 影响 。 有 风 障 工 况 上 游 边界 不 同 、 下 游 边界 相同 地 形 的 比较 (地 形 1c, 2c 和 3c)， 风 
障 在 天 线 位 置 前 弱 风 速 的 大 小 同样 与 进入 核心 区 域 的 风 攻 角 呈 正 相 关 。 上 游 边 界 相同 、 下 游 边界 不 同 
的 3 组 地 形 ( 地 形 3c, 3b 和 3a) ， 下 游 边 界 最 高 点 高 程 差 很 小 ， 对 天 线 位 置 削弱 风速 的 影响 差异 也 很 
c 小 ， 调 控 后 天 线 位 置 的 风速 分 别 为 7. 82 m/s, 7. 81 m/s 和 7. 78 m/s, 
rz ZEA S. AL 10 以 及 表 2 分 析 南 方向 来 风 边 界 地 形 对 风 场 调控 的 影响 。 有 风 障 工 况 上 游 边界 不 
QO 同 、 下 游 边 界 相 同 地 形 的 比较 (地 形 3c, 3b 和 3a) ， 由 于 没有 和 斜坡 地 形 的 影响 ， 风 场 调控 后 天 线 位 置 
处 的 风速 表现 出 与 人 流风 攻 角 的 相关 性 。3c, 3b 和 3a 地 形 的 入 流风 攻 角 分 别 为 1.62^, 2. 18^ 和 
1. 43”， 天 线 位 置 的 风速 分 别 为 7.20 m/s, 7. 05 m/s 和 7. 96 m/s， 风 攻 角 越 大 ， 调控 后 的 风速 越 小 。 
有 风 障 工 况 上 游 边 界 相 同 、 下 游 边界 不 同 地 形 的 比较 (地 形 1e, 2c 和 3c)， 由 于 下 游 没 有 回流 旋涡 ， 
下 游 边界 地 形 对 天 线 位 置 风 场 调控 的 影响 主要 表现 在 地 形 对 流速 的 阻碍 上 ， 下游 边界 最 高 点 海拔 越 
高 ， 风 场 调 控 后 天 线 位 置 的 风速 越 小 。 

2.4 讨 论 

对 于 目标 区 域 相对 水 平 的 地 形 工 况 ， 风 障 遮 蔽 区 降低 风速 的 大 小 主要 受 上 游人 流风 攻 角 的 影响 ， 
入 流风 攻 角 越 大 ， 风 障 遮 蔽 区 降低 的 风速 值 武大 。 另 外 ， 下 游 地 形 海拔 越 高 ， 对 风流 的 阻碍 越 大 ， 风 
障 遮 项 区 降低 风速 的 值 也 越 大 。 如 果 下 游 地 形 高 度 比 风 障 高 度 低 ， 则 下 游 地 形 对 风速 折 减 系数 的 影响 
不 明显 ， 如 图 8。 

台 址 的 斜坡 地 形 增 加 了 流 场 的 复杂 程度 ， 和 斜坡 南 高 北 低 ， 增 强 了 南边 边界 地 形 对 风速 折 减 系数 的 
影响 ， 同 时 削弱 了 北边 边界 地 形 对 风速 折 减 系数 的 影响 ， 如 图 4 和 图 8。 上 游 边 界 地 形 的 海拔 越 高 ， 
对 下 游 区 域 的 挡 风能 力 越 强 ， 但 下 游 树 立 的 风 障 ， 其 减 风能 力 一 定 程度 上 有 所 削弱 ， 且 风速 折 减 系数 
与 地 形 高 度 非 线 性 相关 。 因 此 ， 在 进行 风 障 风 场 仿真 时 ， 尤 其 是 相对 复杂 的 地 形 工 况 ， 上 游 边界 要 尽 
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量 延 伸 到 相对 较 高 山体 的 外 围 ， 以 保障 风 场 仿真 的 可 靠 性 。 表 1 中 北方 向 来 风 ，1a 地 形 在 天 线 位 置 
的 风速 折 减 系数 为 40. 65 96, 2a 地 形 在 天 线 位 置 的 风速 折 减 系数 为 47. 31%， 差 值 大 于 6%。 
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图 8 台 址 水 平地 形 顺风 向 风速 折 减 系数 分 布 
Fig.8 Distribution of reducing wind speed coefficients downwind on the horizontal terrain of the site 
R2 人 台 址 水 平地 形 天 线 位 置 风速 折 减 系数 
Table 2 Ranking of reducing wind speed coefficients for antenna location on horizontal terrain of site 
Wind from north Wind form south 
Wind speed Wind speed Reducing Wind speed Wind speed Reducing 
Terrain without wind with wind wind speed Terrain without wind with wind wind speed 
number barrier barrier coefficient number barrier barrier coefficient 
/ (m/s) / (m/s) / (96) / (m/s) / (m/s) / (96) 
3a 22.40 7.78 65. 26 1b 21. 99 6.75 69. 31 
3b 22. 46 7.81 65. 21 le 22. 04 6. 90 68. 68 
3c 22. 46 7. 82 65. 19 2b 22. 23 7.01 68. 45 
2a 22.25 8.21 63. 09 3b 22.31 7. 05 68. 39 
2b 22.31 8.24 63. 05 2c 22.26 7.17 67. 80 
2c 22.31 8. 25 63. 03 3c 22. 33 7.20 67. 73 
la 22. 09 9. 96 54. 93 la 21.18 7. 65 63. 88 
1b 22. 16 9. 99 54. 92 2a 21.45 7.92 63. 07 
le 22. 16 9.99 54. 90 3a 21.54 7.96 63. 05 
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Fig.9 Downwind wind flow streamline on the horizontal terrain of the site with northward wind 
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B10 南方 向 来 风 台 址 水 平地 形 顺风 向 风 场 流 线 


Fig. 10 Downwind wind flow streamline on the horizontal terrain of the site with southward wind 
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分 析 图 4 可 知 ， 台 址 实际 斜坡 地 形 ， 双 风 障 在 北方 向 来 风 的 挡 风 效率 要 远 弱 于 南方 向 来 风 的 挡 风 
效率 。 因 此 ， 风 障 的 实际 布置 ， 北 边 布置 的 风 障 高 度 要 高 于 南边 风 障 的 高 度 。 比 较 图 6 和 图 10， 相 
同 孔 隙 率 的 风 障 ， 和 斜坡 地 形 天 线 位 置 的 回流 旋涡 明显 ， 而 水 平地 形 天 线 位 置 没 有 回流 。 触 发 流动 分 离 
现象 产生 的 条 件 ， 除 了 地 形 因素 ， 还 与 风 障 的 高 度 、 孔 隙 率 等 参量 有 关 。 在 确定 的 地 形 工 况 条 件 下 ， 
可 以 通过 改变 风 障 高 度 、 孔 隙 率 来 调节 风 障 下 游 的 流 场 。 


3 结 论 


本 文 以 奇 台 射电 望远镜 台 址 地 形 为 研究 对 象 ， 利 用 数值 模拟 方法 研究 了 不 同 边界 地 形 对 台 址 风 场 
调控 的 影响 。 为 了 比较 台 址 斜坡 地 形 对 风 场 流动 的 主导 影响 以 及 明确 边界 地 形 对 风 场 调控 的 影响 关 
系 ， 我 们 构建 了 台 址 水 平地 形 对 照 组 。 对 不 同 地 形 仿 真 结果 的 分 析 表 明 ， 在 无 风 障 工 况 中 ， 无 论 台 址 
是 斜坡 地 形 还 是 构建 的 水 平地 形 ， 天 线 位 置 的 风速 大 小 主要 与 上 下 游 边界 地 形 的 最 高 海拔 有 关 ， 海拔 
越 高 天 线 位 置 的 风速 越 小 。 在 有 风 障 工 况 中 ， 目 标 区 域 为 水 平地 形 时 ， 风 障 遮 蔽 区 降低 风速 的 大 小 主 
要 与 上 游 入 流风 攻 角 有 关 ， 入 流风 攻 角 越 大 (无 风 障 工 况 下 )， 风 障 庶 项 区 降低 的 风速 值 越 大 。 而 风 
障 下 游 ， 地 形 海拔 越 高 ， 风 障 遮 珊 区 降低 风速 的 值 也 越 大 。 台 址 的 斜坡 地 形 增加 了 流 场 的 复杂 程度 ， 
在 进行 风 障 风 场 仿真 时 ， 上 游 边界 要 尽量 延伸 到 相对 较 高 山体 的 外 围 ， 以 保障 风 场 仿真 的 可 靠 性 。 当 
风 障 布设 的 位 置 确定 时 ， 可 以 通过 改变 风 障 高 度 、 孔 际 率 等 参量 来 调节 风 障 下 游 的 流 场 。 
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Research on Wind Flow Regulation Influenced by the Mountain Terrain 
at Large Aperture and High Precision Radio Telescope Site 
Based on Numerical Simulation 
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Abstract: Optimizing the wind environment at the site is important to improve the effective observation 
time of large-aperture and high-precision radio telescopes in the middle and high frequency bands. Through 
wind barrier regulating wind flow that can effectively reduce the impact of wind load on the telescope, the 
layout design of wind barriers is not only related to the height, porosity and other parameters of wind barriers , 
but also consider the impact of the terrain on wind shield. However, there are few researches on this area. 
Taking the QiTai radio Telescope ( QTT) site as the research object, the influence of different terrains on 
regulating wind flow is carried out based on numerical simulation. The simulation result shows that the wind 
speed at the telescope location without wind barriers is mainly related to the highest elevation of the upstream 
and downstream boundary terrain. Under the condition with wind barriers, if the telescope is located in 
horizontal terrain, the magnitude of reducing wind speed in the shelter area is related to the wind attack angle 
of the upstream incoming wind and the elevation of the downstream terrain. The actual slope terrain of the site 
increases the complexity of the wind flow, and when doing wind flow simulation for wind barriers, the 
upstream boundary should be extended to the outside of the relatively high mountain as far as possible. The 
simulation result of different boundary terrains show that the difference of reducing wind speed coefficient at 
telescope position reaches a maximum of 6%. This study can provide a reliable reference for terrain modeling of 
regulating wind flow simulation. 


Key words: regulating wind flow; radio telescope site; wind barrier; terrain; numerical simulation 


